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Abstract 

This paper studies a spatial queueing system on a circle, polled 
at random locations by a myopic server that can only observe cus- 
tomers in a bounded neighborhood. The server operates according 
to a greedy policy, always serving the nearest customer in its neigh- 
borhood, and leaving the system unchanged at polling instants where 
the neighborhood is empty. This system is modeled as a measure- 
valued random process, which is shown to be positive recurrent under 
a natural stability condition that does not depend on the server's scan 
radius. When the interpolling times are light-tailed, the stable system 
is shown to be geometrically ergodic. The steady-state behavior of the 
system is briefly discussed using numerical simulations and a heuristic 
light-traffic approximation. 
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1 Introduction 

This paper studies a spatial queueing system where customers arrive to 
random locations in space that are a priori unknown to the server. The 
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server operates sequentially in time by scanning a bounded neighborhood of 
a randomly chosen location, and serving the nearest customer it observes. 
After a service completion or an unsuccessful scan, a new scan is performed. 
Such models are motivated by applications, such as disk storage or wireless 
sensor networks, where the time required for the server to find a customer 
may have a major impact on the customer's sojourn time in the system. 

To simplify the analysis, we will restrict to systems where service times 
are negligible compared to scanning times, by assuming all service times 
to be identically zero. Observe that in the opposite case where the service 
times dominate, the number of customers in the system behaves like the 
standard single-server queue. A further simplification is to assume that the 
customers arriving during a scanning period are ignored from the ongoing 
scan. We expect that this assumption is valid when the scan radius is 
small. The resulting system can be rephrased as follows: A server polls 
the system sequentially at random locations, and either serves the nearest 
observed customer immediately, or leaves the system unchanged if there are 
no customers within the scan radius. Instead of a single polling server, we 
may alternatively think that there is an infinite stream of servers, each of 
which can only perform one scan during its operational lifetime. 

The service policy is greedy in the sense that the server always aims 
to serve the nearest customer. We call this model greedy polling server, 
to distinguish it from spatial queueing systems where the server travels 
in space towards a nearest customer (greedy traveling server). Although 
a natural choice for many applications, the greedy policy is hard to an- 
alyze mathematically, as confirmed by earlier studies on greedy traveling 
servers [2, 3, 6, 12, 17, 18]. Coffman and Gilbert [(>] conjectured that 
the greedy traveling server on a circle is stable when the traffic intensity 
is less than one, regardless of the server's traveling speed. Kroese and 
Schmidt [17] proved this statement for several nongreedy policies, and Foss 
and Last [12, 13] proved it for greedy traveling servers on a finite space (see 
also Meester and Quant [21]). Further, Altman and Levy [2] analyzed the 
stability of the so-called gated-greedy policy on convex spaces, and Altman 
and Foss [1] derived stability conditions for nongreedy randomly traveling 
servers on general spaces. Despite these affirmative results for closely related 
systems, the stability of the greedy traveling server on a circle still remains 
an open problem. 

The distribution of customers in continuum spatial queueing systems 
with a nongreedy server is now relatively well understood (Kroese and 
Schmidt [17]; Eliazar [n, !)]; Kavitha and Altman ['"]), whereas analyt- 
ical results related to greedy traveling servers are scarce. Coffman and 
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Gilbert [fi] showed that the greedy travehng server on a circle closely re- 
sembles a cyclic policy in heavy traffic (see also Litvak and Adan ['2()] for 
a similar result). Bertsimas and van Ryzin [ ] found two lower bounds for 
the mean sojourn time in the system, which are valid for all travel poli- 
cies. Kroese and Schmidt [18] derived second-order approximations for the 
number of customers and workload in light traffic. 

The analytical challenges related to greedy traveling servers suggest that 
greedy polling servers may be difficult to analyze as well. This is why we set 
a modest research goal for this paper, namely, to characterize the stability 
of the system. This problem may be viewed as determining the system's 
throughput capacity, defined as the maximal sustainable arrival rate for 
which the number of customers in the system remains stable (stochastically 
bounded). Foss [11] recently presented an open problem, attributing it to 
V. Anantharam, conjecturing that the greedy polling server on a circle is 
stable if the arrival rate is less than the polling rate, regardless of the server's 
scan radius. 

In this paper we prove Anantharam's conjecture [11, Section 3.2] by 
showing that the greedy polling server on a circle is stable if and only if 
the polling rate exceeds the arrival rate. The proof is based on presenting 
the system as a measure- valued Markov process, and developing a novel 
quadratic Lyapunov functional on the space of finite counting measures for 
which the measure- valued process has negative mean drift for large customer 
configurations. Besides incrementing the collection of known provable facts 
on spatial queueing systems with greedy service, our analytical results may 
be interesting in other application areas. For instance, the greedy polling 
server can be viewed as a spatial birth-and-death process. Spatial birth-and- 
death processes have usually been studied in the case where all individuals 
have a constant death rate, see for example Ferrari, Fernandez, and Gar- 
cia [10]; and Garcia and Kurtz [14], who give sufficient conditions for sta- 
bility in terms of the birth rates. The greedy polling server differs from the 
above birth-and-death processes in that the death rates of individuals are 
governed by the Voronoi tessellation generated by the customer locations. 
Borovkov and Odell ['] have recently studied a class of spatial-temporal 
point processes based on Voronoi cells, where the number of individuals is 
assumed constant over time. We expect that the quadratic Lyapunov func- 
tional presented in this paper may turn out useful in studying the ergodicity 
of more general spatial birth-and-death processes. 

During the final writing stage of this article, we came across an inter- 
esting recent work of Robert [2-)], who considers the same problem from a 
different point of view. Using entirely different techniques (a stochastically 
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monotone construction of a stationary solution), he proves a weaker form of 
stability, stating that the system has a limiting distribution for which the 
number of customers is finite almost surely. Our results based on Foster- 
Lyapunov drift criteria allow to prove stronger forms of stability, such as 
positive Harris recurrence and geometric ergodicity, depending on the tail 
behavior of the interpolling time distribution. 

The rest of the paper is organized as follows. In Section 2 we describe the 
system as a measure- valued Markov process and derive formulas for its tran- 
sition operators. Section 3 shows the positive recurrence of the system, and 
Section 4 is devoted to geometric ergodicity. In Section 5 we illustrate the 
system dynamics with numerical simulations complemented with a heuristic 
approximation of the system in light traffic. Section 6 concludes the paper. 

2 System description 
2.1 Notation 

The server operates on the circle S = {{xi,X2) G : + = d"^} of 
circumference i > 0, where the distance d{x, y) between points x and y is 
defined as the length of the shortest arc connecting them. The state space of 
the system is the set M^{S) of finite counting measures ( on S, so that C(-S) 
denotes the number of customers located at -B C S*. The elements of M^{S) 
will also be called configurations, and the zero counting measure is called 
the empty configuration. The total number of customers in configuration 
C is denoted by ||C||v = CiS); this quantity also equals the total variation 
of (. We equip the space M+(S') with the sigma-algebra generated by the 
maps C I— )• CiB), with B being a Borel set in S. For real functions / on S 
we denote 



where x G C is shorthand for C{{x}) > 0- For more background, see for 
instance Daley and Vere- Jones [ ]. 

2.2 Population process in continuous time 

Customers arrive to the circle S at uniformly distributed random locations. 
Assuming that the arrival locations and interarrival times are independent, 
the spatial-temporal arrival process can be described as a Poisson random 
measure on x S with intensity measure \dtm{dx), where A > de- 
notes the mean arrival rate, dt the Lebesgue measure on R_|_, and m{dx) 
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the uniform distribution (Haar measure) on S. The server polls the circle 
sequentially in time by scanning a neighborhood of radius r > of a ran- 
domly chosen location, and either serving the nearest observed customer 
immediately, or leaving the system unchanged if the scanned neighborhood 
is empty. Hence the probability that a customer located at x E C is served 
during a polling event is equal to 

m{BrnT(;{x)), 

where Br{x) denotes the open r-ball centered at x, and 

T^{x) = {yeS: d{x,y) < d{x',y) Vx' G (} 

denotes the Voronoi cell of a point x with respect to configuration (If 
there are many customers located at x, each of them is served with equal 
probability.) Assuming that the interpolling times are independent and 
exponential, say with parameter ^, and independent of the arrival process, 
the customer population in the system is described by a Markovian spatial 
birth-and-death process [ ] with generator 

MC) = A / (/(C + 5.)-/(C))m(dx) 
Js 

+ / (/(C - <5x) - /(C)) m{Brix) n r^(x)) Cidx), 
Js 

where 6x denotes the Dirac measure at x. 

2.3 Population process at polling instants 

The stability regions of many ordinary queueing systems are insensitive to 
the shape of the service time distribution. We will show in Section 3 that 
analogously, the shape of the interpolling distribution does not affect the 
stability of our model. To show this claim, we will from now on assume 
that the interpolling times follow a general distribution G on M_|_ with a 
finite mean. To characterize stability in terms of finite mean hitting times 
into small sets (see Lemma 2.1 in Section 2.4), it is sufficient to study the 
discrete-time population process W obtained by sampling the system at 
polling instants, so that Wt{B) denotes the number of customers in B C S 
just after the t-th polling instant, t G Z_|_ (we assume that also the initial 
state Wq is observed just after a polling instant). Observe that if the mean 
hitting time of W into a set is finite, then the same is true also for the 
continuous-time population process, because the mean interpolling time is 
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assumed finite. Therefore, we will from now on only analyze the discrete- 
time population process W. 

The population process is a discrete-time Markov process in M_|-(S'). 
Given an initial state C, G M_|_ (S) , we denote the one-step transition operator 
of by 

AfiO = f{Wi) = E (/(Wi) I W-o = , 

where / is a bounded or positive measurable function on M+(S'). The 
associated probability kernel will be denoted by 

P{C,B) = AlB{0 = PdWieB), BcM+iS). 

Because the polling locations are independent of the arrivals, we can 
decompose the transition operator according to 

A = A^o Ap, 

where the arrival operator A^, and the polling operator Ap are defined as 
follows. The operator A^ acts on bounded or positive measurable functions 

by 

^a/(C) = E^o/(C)GA(n), (1) 

n>0 

where 

^o/(C)= / f{C + S,)m{dx) 
Js 

corresponds to adding one new customer to a uniform random location, and 

is the probability that n customers arrive during an interpoUing time. The 
operator Ap is defined by 

^p/(0 = f{x){l - kriO) + E - S.) m{Br{x) n r^(x)), (2) 

where 

kriO = m{Ua,ecMx)) (3) 

is the probability that the server finds a customer during a scan targeted 
into configuration 
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2.4 Irreducibility and aperiodicity 

The following result shows that the Markov process W describing the cus- 
tomer population in the system may become empty with probabilities uni- 
formly bounded away from zero. As a consequence, the process satisfies the 
irreducibility and aperiodicity properties summarized below (see Meyn and 
Tweedie ['^'^] for details). 

Lemma 2.1. The Markov process W is (p-irreducible and strongly aperiodic, 
where (j) is the Dirac measure on Mj^{S) assigning unit mass to the zero 
counting measure on S. Moreover, the level sets of the form Cn = {C £ 
M+(5) : llCllv <n],n> 0, are small for W. 

Proof. Consider an initial configuration Q with ||C||v = k < n customers. 
Then the probability that the system is empty after n polling instants is 
greater than the probability that no customers arrive during the first n 
interpolling times and the server finds a customer at each of the k first 
polling events. Because the probability of finding a customer in a nonempty 
system is at least m{Br), it follows that 

P"(C,{Co})>e"m(i?,)'=>e"m(i?,,r, 

where Co denotes the empty configuration and e = f^^ e~^^ G{ds) is the 
probability that no customers arrive during an interpolling time. Denoting 
/i = e^m{Br-)^(j), we may thus conclude that 

inf P^{C,B)>^l{B) 

for all measurable B C M+(S'). Thus the level set C„ is small [22, Section 
5.2]. Further, the inequality ^((^OjICo}) ^ e implies that W is strongly 
aperiodic [22, Section 5.4]. □ □ 

3 Positive recurrence 

This section is devoted to deriving the main stability results (Theorems 3.8 
and 3.11), which together imply that the system is positive recurrent if and 
only if the arrival rate of customers is strictly less than the polling rate, 
regardless of the server's scan radius. We start in Section 3.1 by discussing 
why it is not sufficient to analyze the mean drift of the population size, 
and introduce in Section 3.2 a quadratic functional for which the mean 
drift analysis works. Section 3.3 discusses a key interpolation inequality 
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that is applied in Section 3.4 to show that the mean drift with respect to 
the quadratic functional is negative for large configurations. Section 3.5 
summarizes the behavior of the unstable system. 

3.1 Mean drift with respect to population size 

A common method to prove the stochastic stability of a queueing system 
is show that the mean drift of the system with respect to the number of 
customers is strictly negative for large configurations. To see why this ap- 
proach is not sufficient for the greedy polling system in this paper, denote 
the number of customers in configuration by = ||C||v) and recall that 
the mean drift of the system with respect to h is defined by 

Dh{Q) = Ah{0-h{0, 

where A is the one-step transition operator of the system. Recall the de- 
composition of ^4 = o in Section 2.3, and observe that 

^a/i(C) = MC)+AS1, 

where si = J s G{ds) is the mean interpolling time, and 

A^h{0=h{C)-kr{0, 

where k^- is the probability of a successful scan defined by (3). As a conse- 
quence, 

Dh{0 = \si-A^kr{0. 

Consider a configuration C, = n5x, where n customers are located in a 
single point x £ S. Then kr{C) = m{Br), so by conditioning on whether 
customers arrive or not during a polling instant, we find that 

A.kriC) < A:,(C)Ga(0) + (1 - Ga(0)) 
= l-GA(0)(l-m(i?,,)), 

where Ga(0) = / e-^'Gids). Hence 

Dh{C) > Asi - 1 + Ga(0)(1 - m{Br)). 

Because the right side above does not depend on n, we see that Dh{C,) can 
be strictly positive for arbitrarily large configurations, if Asi > 1 — Ga(0)(1 — 
m{Br))- Hence the mean drift with respect population size cannot be used 
to show that the system is stable whenever Asi < 1. 
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3.2 Quadratic energy functional 

Let M{S) be the space of signed counting measures on S, that is, measures 

of the form ( = X^ILi ^i^Xi, where Zj G Z and Xj G S. The space M{S) is 

a normed vector space with the total variation norm ||C||v = J2i 

the subspace of finite positive counting measures M^{S) is a convex cone in 

M{S). 

Given < a < £/2, define for C, r/ G M(5), 

(C,?/)a= / / {a- d{x,y))+C{dx)r]{dy), 
Js Js 

and denote 

iiciia = \/(c:c):. 

Lemma 3.1. The map {Cv) (C;^)a symmetric, bilinear, and positive 
semidefinite. The map (" i— >• ||C||a a seminorm on M{S) satisfying 

IICl|a<||Cl|v (4) 

for all C G M(S'). 

Proof. Clearly, (C)^)a is bilinear and symmetric, so we only need to show 
positive semidefiniteness. We shall prove this using a probabilistic argu- 
ment, representing the bilinear functional as an expectation of a random 
quantity. Observe that (a — d{x,y))+ = 771(5^/2(2^) -^0/2(2/))) where m 
denotes the uniform probability distribution on S. Hence by changing the 
order of integration we see that 

(C,??)a= / / m{Ba/2ix)nBa/2iy))Cidx)r]{dy) 
Js Js 

= j j f l{x€Ba/2{zmyGBa/2{z))m{dz)C{dx)r]{dy) 
'J S "J S J s 

= [ [ [ l{z B,/2{xmz € B,/2{y))C{dx)v{dy)m{dz). 
Js Js Js 

Hence 

(C,r/)a = EC(i?./2(f/)Mi?a/2(^^)), 

where C/ is a random variable uniformly distributed on S. Especially, 

(C,C>a = EC(i?„/2(^))^ (5) 

which shows that (CiC)a > for all C £ M{S). As a consequence, ( 1— )• ||C||a 
is a seminorm (see for example Rudin [24, 4.2]). The representation (5) 
further shows the validity of (4). □ □ 
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Remark 3.2. When a is chosen so that i/a is not an integer, one could 
perhaps strengthen the statement of Lemma 3.1 by showing that (C;'?)a is 
an inner product on M{S). However, in the sequel we will only need the 
fact that I Id la is a seminorm. 

Lemma 3.3. For any integer n and any configuration C G M^{S), there 
exists a closed ball B in S with diameter n^^ such that 

aB)>n-'\\C\U. (6) 

Proof. Given an integer n, cover the unit circle with closed balls Bi, . . . , Bn, 
each having diameter n~^. Then 



v< VC(^*)<nmaxC(Bi), 



i=l 



which shows that maxj > n "^||C||v D D 

Lemma 3.4. For any a > for all ( G M^{S), 



l + 2/a 



<||CI|a<^/^||CI|v. (7) 



Proof. The upper bound in (7) follows directly by observing that (a — 
< a for all x and y. To prove the corresponding lower bound, 
let n be an integer such that 2/a < n < 2/a + 1, and use Lemma 3.3 to 
choose a closed ball B with diameter such that (6) holds. Because 
(a — d{x, y))+ > a — n"^ > a/2 for all x,y ^ B, it follows that 

IICI la > ^ J^{a - d{x, y)U Cidx) ({dy) > ^C{B)^ 

Because n < 2/a + 1, the lower bound now follows using (6). □ □ 



3.3 Interpolation inequality 

This section is devoted to proving a key inequality (Lemma 3.7), which 
is needed for analyzing the mean drift of the system with respect to the 
seminorm ||C||a- For a point x on the circle S and a positive real number a, 
we denote by x + a the point on the circle obtained by traveling distance a 
from x anticlockwise on the circle, and x—a the corresponding point obtained 
by traveling in the clockwise direction. Moreover, we denote by [u, v] the 
closed arc formed by drawing a line from u to v moving anticlockwise on the 
circle, and by {u,v) the interior of [u,v]. 
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Lemma 3.5. Let < a < i/2. Then for all x e S and for all C G M+{S) 
having atoms at x — a, x, and x + a, 

Y,{a-d{x,y))+m{T^{y))=a\ (8) 

Remark 3.6. Equation (8) may interpreted in terms nearest-neighbor in- 
terpolation as follows. Given a bounded measurable function / on S, define 
its nearest-neighbor interpolant (with respect to configuration C,) by 

kf{z) = Y.f{y)l{z(^T^{y)). 

Then 

" i^f{z)m{dz) = J2fiy)M^dy))- 



L 



Given ( and x as in Lemma 3.5, define f{z) = (a — d{x,z))^. Then the 
left side in (8) equals Jg I,^f{z) m{dz), and a short calculation shows that 
Jg f{z) m{dz) = a^. Hence (8) can be rewritten as 



Ic_f{z)m{dz) = I f{z)m{dz), 
Js 

stating that the interpolation error made in replacing / by /^/ is zero (see 
Figure 1). 

Lemma 3.5. Assume first that x is the only atom of C, in (x — a,x + a). Then 

m{Tc_{x)) = m{[x — a/2, x + a/2]) = a, 

which implies (8). 

Assume next that (8) holds for some configuration ( having atoms at 
X — a, X, and x + a. We shall show that the same is true also for (' = 
C + 6z, where z G S is arbitrary. We only need to consider the case where 
z E {x — a,x + a) such that z ^ (, because otherwise the sum on the 
left side of (8) remains unaffected. Assume without loss of generality that 
z £ {x — a,x), the other case being symmetric, and let 2/1,2/2 S [x — a, x] be 
the points of C lying nearest to z in the anticlockwise and clockwise direction, 
respectively. Then 

m{T(;'{z)) = ^{d{yi,z) + d{z, ^2)), 

"i(rc'(yi)) = "^(rc(2/i)) - ^d{z,y2), 

m(r^/(2/2)) = m(rf(y2)) - ;^d(yi,2;). 
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a - d{y,x) 



X — a 



rc(y) 



X 



X + a 



Figure 1: Nearest-neighbor interpolation of 2 i— )• (o — d(x, 2))+. 

Denote by g{x,Q the left side of (8). Because the other atoms of C 
remain unchanged, it follows that 



9{x,C') -g{x,0 = ^{a-d{z,x)){d{yi,z) + d{z,y2)) 

- ^(« - d{yi,x))d{z, 2/2) - d{y2,x))d{yi,z). 

Because d{yi,x) = d{yi,z) + d{z,x) and d{y2,x) = d{z,x) — d{y2,z), the 
terms on the right side of the above equation cancel out each other, so we 
may conclude that g{x,(') = g{x-,C)- Hence (8) holds also for C' = C + ^z, 
and the proof is complete by induction. □ □ 

Lemma 3.7. Assume that < a < min(-^/2, 2r). Then for all C G 
and for all x ^ C,, 



^(a - d{x, y))+m{Br{y) D T^iy)) > 



(9) 



Proof. Denote the left side of (9) by g{x, C,). Because a < |, the sum on the 
left of (9) runs over the locations y S C such that y E [x — a,x + a]. 
Let C' = C + ^x~a + ^x+a- We shall first show that 



g{x,C)>g{x,C). 



(10) 
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Denote by yi and yr the atoms of C in [x — a,x + a] located nearest to 
X — a and x + a, respectively. Then (see Figure 2) r^/(y) = r^(y) for all 
y G C 1^ (yhUr), and moreover, 

r^'(y/) C rf(yO and rf/(yr) C T(^{yr). 

Because the points y = x — a and y = x + a contribute nothing to the sum 
on the left side of (9), we may conclude that (10) holds. 

Further, because a < 2r, it follows that r^'(y) n Br{y) = r^/(y) for all 
y G C' 1^ (^ ~ ^) 2; + a). Hence Lemma 3.5 shows that 

9{x, C) = ^{a- d{x, y))+m(r^(y)) = a^. 
In light of (10), this implies the claim. □ □ 



x+'a 



X — a 



yi 



Figure 2: Voronoi cells of yi and yr- 



3.4 Positive Harris recurrence 

The following result shows that the system is stable under the natural condi- 
tion Asi < 1. Recall that si and S2 denote the first and the second moments 
of the interpolling time distribution, respectively. For a converse to Theo- 
rem 3.8, see Theorem 3.11 in Section 3.5. 

Theorem 3.8. Assume Xsi < 1 and S2 < 00. Then for all values of the 
scan radius r > and the circle circumference i > 0, the population process 
W is positive Harris recurrent with stationary distribution vr such that 



I ||C||vVr(dC)<oo. 



Moreover, for all initial states ( G M^{S), 



sup 

9:|9|<1 + I 



gdTT 



as t 



00. 
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Remark 3.9. Meyn and Tweedie call such a process /-ergodic with /(C) 
1 + 



The proof of Theorem 3.8 is based on the following Foster-Lyapunov 
bound on the mean drift of the system with respect to the functional h{Q = 

(COa. 

Lemma 3.10. Assume that si, S2 < oo, and letO<a< min(^/2, 2r). Then 
the mean drift of the system with respect to h{Q = (C;C)a satisfies 

Dh{C) < -Cl||C||v + C2 

for all Q G Mj^{S), where 

ci = 20^(1 - Asi), 

C2 = a(l + Asi) + a^{\^S2 - 2\si). 

Proof. RecaU that Dh{C,) = Ah{C) - h{C), where A = A;, o Ap, and the 
transition operators A^, and Ap are defined by (1) and (2). Observe that the 
transition operator corresponding to n arrivals satisfies 



1=1 

n n n 

= hiO + 2 E(C, 5x,). + 6x, , ^ 6x,). 

i=l i=l j=l 

= h{C) + 2n E(C, 5xi >a + n E{Sx, , dx^ >a + (n^ - n) E{5x, , 6x2 >a, 

where Xi are independent and uniformly distributed on S. Because the 
uniform distribution m on S" is shift-invariant, 

(a - d{x,y))+m{dy) = 

JS 

for all X G S, so that 

E(C,%>a= / y2{a-dix,y))+C{{y})midx)=a' 



and 



E{Sxi,Sx2)a.= / / {a - d{x,y))+m{dx)m{dy) = a"^ 
Js Js 
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Hence 

= h{C) + 2no2||C| Iv + na + (n^ - n)a2. 

Because 

•^'*+ n>0 

we find using (1) that 

AMO = HO + 2Asia2||C||v + Asia + A^ssa^, 

where Sk = f G{ds) denotes the k-th. moment of the interpolling distri- 
bution. 

To calculate Aph, observe that 

/i(C-<5,)-MC) = -2(C,5..)a + a, 

which shows that 

AphiC) = HO + akriO - 2 J^(C, 5.,)^m{Br{x) n T(;{x)), 

where kr{C) is the probability of a successful scan, defined by (3). Lemma 3.7 
implies that for all ( G M+(S') and all x £ 

^(a - d{x, y))+ m{Br{y) n T^iy)) > a^, 

I/6C 

SO especially, 

J^(C, Sy), m{Br{y) n r^{y)) > a^\\C\\.. 

Hence 

Aph{0<h{C) + a-2a^\\C\\^. 

Because 

Aa/i(C) = HO + 2Asia2||C||v + Asia + A^saa^ 

and 

^a(||-||v)(C) = IICI|v + Asi, 

we find that 

AHC) < HO + 2Asia2||C||v + Asia + A^saa^ + a - 2a2(Asi + ||C||v). 

□ □ 
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Theorem 3.8. By Lemma 3.10, the mean drift of the system with respect to 
v{Q = WQWl for < a < min(£/2,2r) satisfies 

Dv{Q < -Ci||C||v + C2, 

for ci > and C2. Define /(C) = 1 + c||C||v where c = ci/2, and let n be an 
integer such that n > (1 + C2)/(ci/2). Then 

^^;(C) < -/(C) 

for all C G ^+{S) such that ||C||v > n. Moreover, 

Dv{C) + f{0<l + C2 

for all C S M+(5). Lemma 2.1 shows that the system is c/i-irreducible and 
aperiodic, and the level set Cn = {C, ^ M^{S) : ||C||v < n} is small and 
thus petite. Hence the /-norm ergodic theorem of Meyn and Tweedie [22, 
Theorem 14.0.1] shows the claim. □ □ 

3.5 Instability 

Theorem 3.11. If Xsi > 1, then the system is not positive recurrent, and if 
Xsi > 1, then \ \Wt\\v oo almost surely as t ^ oo, regardless of the initial 
state. 

Proof. Denote the population process of the system by W, and let W' be 
the population process of a modified system with scan radius r' = i/2, 
so that the server finds a customer at polling events whenever the system 
is nonempty. Then ||VK'||v equals the number of customers in a standard 
single-server queue with rate-A Poisson arrivals and service times distributed 
according to G, observed just after service completions. It is not hard to see 
that there exists a coupling of W and W so that ||Wj||v > ||W^tllv for ah t G 
Z_|_ almost surely, whenever ||Wo||v > ll^ollv (see for instance Leskela [i'.). 
Theorem 4.8]). When Asi > 1, it is well-known that the mean return time 
to zero of ||VF'||v is infinite [ ]. The coupling then implies that the same is 
true for the process ||Ty||v) which shows that W' is not positive recurrent. 
The second claim follows using the same coupling, because ||VF/||v — >• oo 
almost surely when Asi > 1. □ □ 

4 Geometric ergodicity 

The standard single-server M/G/1 queue is known to be geometrically er- 
godic, when the tail of the service time distribution is light enough (Spieksma 
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and Tweedie [ ])• An analogous result is true for the population process 
W, as the next result shows. 

Theorem 4.1. Assume that Xsi < 1 and the interpoUing time distribution 
satisfies f e^'^G{ds) < oo for some 6 > 0. Then the system is geometrically 
ergodic in the sense that there exist constants q > 0, /3 > 0, and c < oo such 
that 

oo ^ 

<ce^llClla 



i=0 9:|9l<ePll-ll=' 



E^g(Wt)- j gdn 



for all initial states C £ Mj^{S). Moreover, the stationary number of cus- 
tomers is light-tailed in the sense that 



s^WiW- TT{dC) < oo (12) 



for some 7 > 0. 



Before proceeding with the proof of Theorem 4.1, we will show that the 
seminorm ||C||a satisfies a similar Foster-Lyapunov bound (Lemma 4.3) as 
the function (C,C)a = llClla ™ Lemma 3.10. Using the Foster-Lyapunov 
bound for ||C||a, we will then proceed along similar lines as in the proof 
of [22, Theorem 16.3.1] (see also Borovkov and Hordijk [ :]) to bound the 
mean drift of the system with respect to the function e'^"'''!'' for some (3 > 
(Lemma 4.4), which is key to proving Theorem 4.1. 

Lemma 4.2. The function x 1— )• \/l — x satisfies 

\/l — X = 1 — -X + R{x), 
where \R{x)\ < 2-3/2 allxe[-^,^]. 

Proof. Taylor's first order approximation shows that for all x G [—5,5]; 
there exists s G [0, 1] such that 

R(x) = -(l-sx)-3/2x2. 
8 

Because (1 — sx)-^/^ < 2^/^ for all |2;| < ^ and s G [0, 1], the claim follows. 

□ □ 

Lemma 4.3. Assume that Xsi < 1 and S2 < 00, and letO < a < min(£/2, 2r) . 
Then there exist a > 0, b > 0, and an integer n such that the mean drift of 
the system with respect to the seminorm v{C) = ||Cl|a satisfies 

Dv{C) < -a (13) 
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for all C G M+{S) such that ||C||v > n-, and 

Dv{C) < b (14) 
for all C £ M+{S) such that ||C||v < n. 

Proof. Jensen's inequality shows that Av < {Av'^y^'^ holds pointwise on 
M+(S'), so the mean drift with respect to v is bounded by 

Dv < (^t;2)l/2 _ ^ = (^2 ^ L»7;2)l/2 _ ^_ (15) 

Because Asi < 1, we see using Lemma 3.10 that there exist c > such that 
Dv'^{() < — c||C||v whenever ||C||v is large enough. Because ||C||a < IICIIv 
(Lemma 3.1), we see that 

Dv < (f^ - cv)^/^ -v = v(^{l- cv"^)^/^ - 1) 

for all llCllv large enough. Lemma 3.4 shows that v{C) — 00 as ||C||v — ^ co in 
M+{S). Thus, for all C G M+{S) such that ||C||v is large enough, cv'^ < \, 
and using Lemma 4.2 we see that 



Dv<v + R{cv'^)^ 



This shows the validity of (13) for a suitable chosen n. 

Lemma 3.10 also shows that Dv^ < ci for all Q E M+(5). Because 
i'(C) ^ llCllv (Lemma 3.1), inequality (15) shows that (14) holds with h = 
(n2+c2)i/2. □ □ 

Lemma 4.4. Assume that Xsi < 1 and f^^ e^* G{ds) < 00 for some 9 > 0, 
and let < a < min(^/2,2r). Then there exist q > 0, /3 > 0, 5 > 0, 
and an integer n such that the mean drift of the system with respect to 
^^/3(C) = exp(/3||C||a) satisfies 

DvpiO < -avpiO (16) 
for all Q € M+(S') such that ||C||v > and 

Dvp{C) + avp{C)<b (17) 
for all C G M+{S) such that ||C||v < n. 
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Proof. Define f/3(C) = e'^"'''!'" for some /3 > 0, to be chosen later. Then the 
mean drift with respect to equals 

Z)t;^(C) = z;^(C) Ec {e'^("^^ll^-ll^»IW - 1} . 

Using a first order Taylor series approximation we see that 

e^^ = 1 + I3t + R{t), 

where the error term is bounded by \R{t)\ < ^^'^t'^e'^^^^ for all t G M. Because 
< s~^e*l*l for all t and aU s > 0, it follows by setting s = that 

This bound implies that 

VfsiCr'DvisiC) < f3Dv{C) + f3^/' e^^^^'^'^^ l|H^il|a-||H^o)l|a|^ 



where we denote v{C) = 

Let us next bound the exponential term. Because ||C||a is a seminorm 
in the space of signed counting measures on S (Lemma 3.1), the triangle 
inequality shows that 

|||W^l||a-||T^0||a| < \\Wl-Wo\U. 

Let us write 

Wi-Wo = r]a- rjp, 

where r]a is a random counting measure describing the arrivals during an 
interpolling time, and rjp is a random counting measure describing the num- 
ber of served customers during the first polling instant {ijp = if the server 
sees no customers and rjp = 5x for some x G Wq + r]a otherwise) . Because 
1 1 • I |a < 1 1 • I |v and ||r/p||v < 1, we find that 

||H^l-^^0||a< ||r/a||v + l. 

Observe next that 



G{ds). 



+ n=0 

A{e(/5+/3'/')-l)s 
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Choose now /? small enough such that \{e^^^^^^^^ — 1) < ^ and e^~^^^^^ < 2. 
Then 

VfsiCr'DvfsiO < PDv{0 + 2/34/3(5(0), (18) 

where G{6) = f e^^G{ds). By Lemma 4.3, Dv{C,) is strictly negative for 
\\C\\v large enough. Hence there exist a > 0,/? > and n such that (16) 
holds for llCllv > n. 

Inequality (18) further shows that 

DvpiC) + avpiO < (f3Dv{C) + 2p^l^G{e) + a) vpiQ. 



for all C, € M^{S). Because ||C||a ^ IICIIv by Lemma 3.1, we have the bound 
vpiC) ^ fi'^" for llCllv ^ n. Inequality (14) in Lemma 4.3 thus shows that (17) 
holds for some h large enough. □ □ 

Theorem 4-1- By Lemma 2.1, the system is (/>-irreducible and aperiodic. Fix 
< a < min(£/2,2r), and choose a > and /3 > as in Lemma 4.4. By 
Lemma 4.4, the function f/3(C) = exp(/3||C||a) satisfies a geometric drift 
condition for the level set C„ = {C G M^{S) : \\C\\v ^ n} for some n large 
enough. The set C„ is small and thus petite by Lemma 2.1. Hence by the 
geometric ergodic theorem of Meyn and Tweedie [22, Theorem 15.0.1], there 
exists a finite number c such that (11) holds. 
Observe next that (11) implies 



giC) - gd'"' 



< ce 



for ah C G M+(5) and for all g : M+{S) ^ R such that \g\ < exp(/3|| - Ha). 
Hence, in light of Lemma 3.4, inequality (12) follows by applying the above 

inequality to the function h{() = exp(7||C||v), where 7 = ( 1^7^ ) □ 



1+2/a 

□ 



5 Population size in steady state 

Having seen that the population process W is positive recurrent for Asi < 1 
and S2 < 00, it would be interesting to find out an explicit expression for the 
stationary distribution of VF, or at least for the stationary distribution of 
the number of customers in the system. Following Kroese and Schmidt [17], 
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we could try the Laplace functional approach. Denote the stationary dis- 
tribution of W by vr, and denote the Laplace transform of the stationary 
population size distribution by 

L{e) = j e-^ll^ll^7r((iC), ^ > 0. 

Then a straightforward calculation shows that 

m = G(A(1 - e-'^)) I e-''ll^llv(l + (e^ - 1)A;,(C)) ^(C), 

where G denotes the Laplace transform of G, and ^^(C) is defined by (3). 
However, solving L{6) from the above equation appears intractable due to 
the nonlinearity of k^- 

To gain insight on the behavior of the number of customers in the sys- 
tem, we have numerically simulated the system for a choice of parameter 
combinations with exponential interpolling times. Figure 3 displays simu- 
lated paths of the population size for A = 0.1 (left) and A = 0.9 (right), 
where si = 1, r = 0.1, and i = 1. The simulations suggest that the sys- 
tem in light traffic is empty a large proportion of time, whereas even for 
moderately heavy traffic (Asi = 0.9), empty system appears a rare event. 




Figure 3: Population size as a function of time in light traffic (A = 0.1, left) 
and heavy traffic (A = 0.9, right). 

In light traffic, we can heuristically argue as follows. Assuming there is 
one customer present in the system, and no new customers arrive, finding 
the customer requires a geometrically distributed number of polling instants 
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with parameter m{Br). Hence, assuming that no new customers arrive, the 
time required to serve the single customer has mean si/m{Br). Thus we 
might expect that the mean number of customers is roughly similar to the 
stationary probability that a renewal on-off process with mean on-time 1 /A 
and mean off-time si/m{Bj.) is in on-state, so that 

To study the accuracy of the above heuristics. Figure 4 plots the approx- 
imation (19) against numerically simulated values of the stationary mean 
population size for varying scan radius r in a system with A = 0.1, si = 1, 
and £ = 1. The plot suggests that the light-traffic approximation is relatively 
accurate for a wide range of scan radii. 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 



Figure 4: Simulated stationary mean population size (dashed line) and the 
light-traffic approximation (solid line) as a function of the scan radius r. 



6 Conclusions 

This paper studied a spatial queueing system on a circle, polled at random 
locations by a myopic server that can only observe customers in a bounded 
neighborhood. Using a novel quadratic Lyapunov functional of the measure- 
valued population process, we showed that the system is positive recurrent 
under a natural stability condition, and proved the geometric ergodicity of 
the system for light-tailed interpolling times. The behavior of the stationary 
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system was discussed in terms of numerical simulations and a heuristic light- 
traffic approximation. The quadratic Lyapunov functional studied in this 
paper appears a promising tool for the analysis of more general spatial birth- 
and-death processes. 
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